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FOREWORD 


The  Project  MAGNET  program,  which  began  in  1951  and  which  is  now  more  than  40  years  old, 
has  primarily  been  tasked  to  globally  collect  high-level,  vector  aeromagnetic  data  to  support 
magnetic  mapping  for  navigational  purposes.  However,  the  unique  remote-sensing  capabilities  of 
the  Project  MAGNET  aircraft  are  only  fully  realized  in  low-level  applications.  In  the  low-level 
mode  of  operations,  this  aircraft  has  been  tasked  to  perform  mineral  resource  evaluations  in 
southeast  Asia,  to  evaluate  hazards  to  submarine  navigation  in  remote  areas  around  the  world,  and 
even  to  search  for  sunken  ships.  Yet,  having  the  ability  to  collect  high-quality  geophysical  data 
sets  from  unique  data  acquisition  platforms  is  of  little  value  without  also  having  developed 
sophisticated  analytic  tools  to  extract  desired  information  from  the  data. 

This  report,  which  concerns  a  low-level  Project  MAGNET  survey  of  the  northern  Juan  de  Fuca 
and  Explorer  Plates,  a  region  of  intense  scientific  interest  due  to  its  volcanic  and  tectonic  activity 
and  due  to  its  proximity  to  the  west  coasts  of  Canada  and  the  United  States,  exhibits  a  generous 
mixture  of  both  high-quality  survey  data  and  sophisticated  mathematical  and  computer  analysis. 
The  results  of  the  analyses  permit  one  to  peer  several  kilometers  below  the  bathymetric  layer  of 
the  ocean  to  see  hidden  basalt  structures  with  unusual  detail.  Occasionally,  one  of  these  basalt 
structures  will  rise  in  the  form  of  a  seamount  well  above  the  usual  oceanic  sediment  coating.  It 
then  becomes  a  hazard  to  submarine  navigation  and  consequently  becomes  of  interest  and 
concern  to  the  navy.  ; 
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1.  INTRODUCTION 


The  Joan  de  Foca  and  Explorer  plates  have  been  extensively  surveyed  over  a  period  of 
several  decades,  primarily  using  shipbome  scalar  magnetometers,  as  noted  by  Raff  and  Mason 
(1961),  Elvers  et  al.  (1974),  and  more  recently  by  Tivey  and  Johnson  (1990).  The  tectonic 
evolution  of  dus  area  is  summarized  by  Botros  and  Johnson  (1988)  and  by  Karsten  and  Delaney 
(!989),  while  the  presoit  geomoiphological  view  of  the  region  is  summarized  by  Finn  (1990). 
The  most  detailed  studies  of  the  region  have  concentrated  on  the  special  relationship  that  exists 
between  the  Axial  and  Cobb  seamounts  and  the  Juan  de  Fuca  Ridge.  This  relationship  is 
described  from  the  magnetic  viewpoint  by  Tivey  and  Johnson  (1990),  from  the  bathymetric 
viewpoint  via  Sea  Marc  I  data  by  Applegate  (1990),  from  the  chemical  viewpoint  by  Rhodes  et 
al.  (1990),  from  the  gravity  viewpoint  by  Hildebrand  et  al.  (1990),  and  from  the  seismic 
reflection  and  refraction  viewpoints  by  Morton  et  al.  (1987)  and  by  White  and  Clowes  (1990). 
The  subduction  of  the  Juan  de  Fuca  plate  undo*  the  western  continental  margin  of  die  North 
American  Plate  has  been  studied  as  far  east  as  the  Wadiington  and  Or^on  Cascade  Mountain 
Range,  using  tdeseismic  P-wave  tomogr^hic  imaging  techniques,  accounts  of  which  are  given 
by  Rasmussen  md  HunqihrQrs  (1988)  and  by  Harris,  Iyer,  and  Dawson  (1991). 

The  original  purpose  for  the  analysis  being  presented  here  was  to  develop  a  practical  method 
of  identifyii^  hazards  to  underwater  navigation  for  the  U.S.  Navy,  using  magnetic  techniques 
duU  take  advantage  of  die  Navy's  unique  aeromagnetic  survey  resource,  the  Project  MAGNET 
RP>3D  Orion  aircraft  The  Project  MAGNET  program  originated  in  1951.  During  the 
program's  lifetime,  several  aircraft  have  been  used  (Colonan,  1992).  The  current  aircraft  has  the 
unique  cq»bility  of  perftnming  large-scale  vector  aeromagnetic  surveys  in  remote  ocean  areas 
in  rektivdy  short  time  ^muis.  The  Juan  de  Fuca  area  provided  a  suffldendy  geologically 


cmnplicated  mturphology  to  test  the  robustn^s  of  any  method  that  might  be  developed.  The 
mott  prominent  hazard  is  the  seamount.  The  basic  assumption  is  that  most  of  these  seamounts 
are  volcanic  in  origin  and  therefore  are  composed  of  basalt  materials  that  have  strong  magnetic 
properties.  The  Juan  de  Fuca  survey  area  includes  the  Heck  seamount  chain,  the  Juan  de  Fuca 
Ridge,  the  Endeavor  Segment,  and  the  Explorer  Ridge,  all  of  which  are  volcanic  in  origin. 

Because  the  aircraft  provides  vector  magnetic  measurements,  it  is  an  easy  matter  to  produce  a 
rectangular-harmonic  magnetic  potential-field  model  from  survey  data  which  have  been 
collected  at  a  single  altitude.  After  appropriate  data  reduction  these  data  can  be  gridded  and 
subsequently  transformed  from  geodetic  to  rectangular  coordinates,  where  the  data  is  then 
amenable  to  Fast  Fourier  Transform  (FFT),  rectangular-harmonic  modeling  techniques,  which 
will  subsequently  be  described.  Technical  aspects  of  rectangular-harmonic  analysis  are 
discussed  by  Alldredge  (1981,  1982).  However,  our  methods  and  procedures  are  somewhat 
different  than  his. 

The  model  presented  here  is  derived  from  survey  data  collected  during  a  detailed,  low-level 
(SOO-ft),  U.S.  Navy  Project  MAGNET,  vector-aeromagnetic  survey  performed  during  the 
months  of  August  and  September  1981.  This  survey  covers  the  ocean  area  bounded  from  4/  N 
to  51^  N  and  from  124**  W  to  130°  W  as  indicated  in  Figure  1.  The  corresponding  bathymetry 
fat  this  region,  taken  from  the  5-minute  gridded  data  base  DBDBS,  is  shown  in  Figure  2.  The 
modd,  in  turn,  is  used  to  recompute  the  magnetic  vector-component  grids  from  which  the  model 
is  derived  and  to  <x>mpute  as  well,  the  nu^etic  gradient-tensor  component  grids  of  the  region. 
These  magnetic  vector  and  gradient-tensor  components  are  combined  via  an  inverse-modeling 
technique  vriuch  enq>loys  the  concept  of  uniformly  magnetized  rectangular  prisms  to  estimate 
the  magnetic  depdi-to-source  for  the  entire  surveyed  region.  As  a  by-product  the  relative 
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Figure  1.  Juan  de  Fuca  Aeromagnetic  Survey  Design;  Altitude:  SOO-ft 
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Figure!.  Juan de Fuct DBDB5 Badiymetry;  Units: 


magnetic  morphology  (i.e.,  the  magnetization  distribution)  of  the  region  is  also  estimated.  The 
rectangular-prism,  inverse-modeling  technique  en:q)loyed  in  this  study  has  a  long  history  with 
many  variations.  The  variation  used  here  is  described  in  detail  by  Quinn  and  Shiel  (1993).  This 
technique  does  not  distinguish  between  remnant  magnetization  and  induced  magnetization. 


DATA  COLLECTION  AND  REDUCTION  PROCEDURES 

The  Juan  de  Fuca  survey  was  conducted  during  the  months  of  August  and  September  1981. 
The  actual  dates  flown  are  given  in  Table  1  along  with  the  project  identification  and  the  flight 
numbers.  These  identifications  are  attached  to  the  data  sets  by  the  Naval  Oceanographic  Office, 
which  conducted  these  surveys  in  cooperation  with  the  Navy's  P3  squadron  VXN-8,  located  at 
Patuxent  River  Naval  Air  Station,  Patuxent  River,  Maryland.  Additional  Project  MAGNET 
flights  in  the  Juan  de  Fuca  area  were  conducted  in  1980  and  1983.  These  flights  were  not  used 
in  this  study  but  are  nevertheless  listed  in  Table  1  as  a  matter  of  completeness  for  the  benefit  of 
future  studies. 

2.1  Instrumentation,  Calibration,  Environmental  Corrections,  and  Processing 

The  1981  survey  was  conducted  at  an  altitude  of  500  feet,  as  determined  by  a  precision 
Rosemont  Barometric  Altimeter  used  by  the  pilot.  The  densely  spaced  East-West  survey  tracks 
were  separated  by  3  luni  except  in  a  magnetically  quiet  zone  located  over  the  trench,  adjacent  to 
the  contmoital  shelf,  where  the  spacing  was  increased  to  6  nmi.  Survey  tracks  flown  for  control 
purposes  in  the  North-South  direction  were  separated  by  20  nmi.  The  effective  along-track 
sample  rate  was  0.5  Hertz,  while  the  speed  of  the  aircraft  was  approximately  240  nmi  per  hour. 
Magnetograms  from  a  nearby  geomagnetic  observatory  at  Victoria,  British  Coliunbia,  were  used 


5 


Tabte  1.  Project  MAGNET  Flights  in  the  Northern  Juon  de  Fuca/Explorer  PUte  Regions 


Project# 

FUght# 

Year 

Month 

Day 

Senior  Scientists 

A32-053 

9102 

1980 

April 

23 

Jesus  Anglero, 
diaries  Gunn  & 

A32  -  053 

9103 

1980 

April 

25 

Mike  Anderson 

A32  - 154 

9411 

1981 

August 

3 

Donald  Wilson, 

Jesus  Anglero  & 

A32  -  154 

9412 

1981 

August 

5 

Steve  Feldman 

A32-154 

9366 

1981 

August 

7 

A32  - 154 

9367 

1981 

August 

10 

II 

A32  - 154 

9413 

1981 

August 

12 

II 

A32  - 154 

9414 

1981 

August 

14 

It 

A32  - 154 

9418 

1981 

September 

3 

It 

C32  -  451 

9429 

1983 

November 

9 

Virgil  Bettencourt  & 
Donald  Wilson 

C32-451 

9430 

1983 

November 

11 

•I 

C32-451 

9431 

1983 

November 

16 

ft 
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to  nuntitor  and  temove  short-term  (i.e.,  <  2  months)  temporal  magnetic  variations.  The  most 
prominent  contribution  to  these  short-term  temporal  variations  is  the  Solar  quiet  (Sq)  Daily 
Variation  (DV)  of  the  genmagnetic  field  and  its  corresponding  induction  field.  Canqrbell, 
Schiffinacher  and  Kroehl  (1989)  have  globally  modeled  the  Sq  variations.  Their  model  is  now 
available  from  the  National  Geophysical  Data  Center  which  is  located  in  Boulder,  Colorado. 
These  Sq  magnetic  variations,  the  ionospheric  currents  that  generate  them,  and  the 
corresponding  induced  magnetic  fields  from  the  Earth's  conducting  crust  and  upper  mantle  are 
described  in  some  detail  by  Campbell  and  Schifimacher  (1988a,  1988b)  and  by  Campbell 
(1989a,  1989b). 

Data  from  the  magnetic  observatory  were  also  used  to  calibrate  the  vector  magnetometer  and 
to  compensate  for  aircraft-generated  magnetic  fields.  Typical  Project  MAGNET 
calibration/compensation  procedures,  which  involve  airswings  over  magnetic  observatories,  are 
described  by  Colenoan  (1992).  A  portable  Vector  Magnetic  Ground  Station  (VMGS)  was 
established  at  McChord  Air  Force  Base  near  Tacoma,  Washington,  also  for  the  purpose  of 
eliminating  teoqtoral  variations  fitnn  the  data.  The  typical  11-  to  12-hour  flights  were 
conducted  at  night  to  further  reduce  the  effect  of  temporal  magnetic  variations  generated  by 
ionosphoic  and  magnetospheric  current  systems  and  their  corresponding  crustal  induction 
effects. 

Navigation  and  attitude  were  determined  using  the  Navy  doppler  radar  satellite  navigation 
systcan  and  a  redundant  set  of  three  inertial  navigation/attitude  systems.  The  RMS  navigation 
error  is  estimated  to  be  on  the  order  of  .1  nmi  RMS,  while  the  attitude  errors  for  Yaw,  Pitch,  and 
Roll  are  estimated  to  be  on  the  order  of  I  arc-minute  RMS.  These  data,  like  the  magnetics  data, 
require  extensive  editii^  for  spurious  noise  spikes,  drift  corrections,  and  time  syncing. 
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GcBmenlly,  the  recorded  time  is  accurate  to  witiun  a  few  milliseconds.  The  Juan  de  Fuca  survey 
took  place  prior  to  the  deployment  of  the  Global  Positioning  System  (GPS).  So,  GPS  navigation 
was  not  available. 

Having  reduced  the  data  for  each  flight,  the  Main  (core-generated)  magnetic  field  was 
removed  using  the  1980  Epoch,  Degree  13  Spherical-Harmonic  Model  (GSFC  12/83)  of  Langel 
nd  Estes  (198S),  which  is  based  on  MAGSAT  data.  The  residual  X,  Y,  and  Z  magnetic  field 
con^nmoits  are  then  essentially  of  crustal  origin.  Data  for  each  magnetic  component  were 
gridded  separately  and  interpolated  as  necessary,  into  lat  x  Ion  cells  of  l.S  arc-minutes  x  1.5 
arc-minutes  (i.e.,  roughly  half  of  the  East- West  track  spacing)  such  that  each  con^nent  consists 
of  256  X  256  grid  points.  The  Project  MAGNET  survey  covered  only  ocean  areas.  It  did  not 
extend  beyond  the  coastline  to  the  adjacent  continental  land  area.  So,  the  extreme  east  and 
noitiieast  portions  of  the  grid  were  zero  filled  in  such  a  way  as  to  allow  interpolation  of  the 
magnetic  residuals  over  the  land  areas  from  the  coastline  to  the  edge  of  the  grid  to  smoothly 
taper  off  to  zero,  thus  avoiding  later  problems  witii  Gibbs  phenomena.  The  gridding  process 
also  has  the  advantage  of  anoothing  out  instrummt  noise,  navigational  errors,  and  other 
environmentally  generated  fields,  not  of  crustal  origin,  that  cannot  otherwise  be  accounted  for 
and  thus  eliminated. 

2Jt  Coordinate  Transformations 

When  the  three  magnetic-corEponent  residual  grids  are  considered  as  a  single  vector  unit, 
each  grid  point  has  associated  with  it  X-,  Y-,  and  Z-componoit  magnetic-field  values 
which  ate  first  transformed  from  geodetic  to  spherical  coordinates  (Cain  et.  al.,  1967)  and 
subsequently  from  spherical  to  rectangular  coordinates  (Alldtedge,  1981).  These  two 
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procedures  involve  not  just  coordinate  transformations  at  the  site  of  each  grid  point,  but  also 
rotations  of  the  residual  magnetic-field  vector  at  each  grid  point. 

The  geodetic  position  (i.e.,  latitude,  longimde,  and  altitude)  of  a  data  point  is  denoted  as 
(X,<pji).  The  same  position  referenced  to  spherical  coordiiutes,  corresponding  to  the  geocentric 
radius,  Ae  colatitude,  and  the  longitude  is  denoted  as  (r,8,(|>).  So,  in  geodetic  coordinates,  the 
residual  ms^etic-field  vector  components,  corresponding  to  X-North,  Y-East,  and  Z-Votically 
Down,  which  forms  a  right-handed  system,  are  denoted  as  B„(X,(p4i),  By(X,(p,h),  and  B/X,(pdt). 
Prior  removal  of  the  Main  field  at  the  mean  survey  epoch  of  1981.2  using  the  GSFC  12/83 
model  and  removal  of  the  short-term  temporal  variations  using  the  observatory  data  eliminates 
the  time  dependence  fiom  the  residual  data  sets.  The  same  residual  magnetic-field 
vector-component  values  referenced  to  spherical  coordinates  have  radial,  colatitudinal,  and 
longitudinal  conqxments  denoted  as  BXr,6,(p),  B,(r,6,9),  and  B^(r,6,(p).  Using  this  notation,  the 
geodetic  to  spherical  coordinate  transformation  used  at  each  grid  point  of  each  of  the  three 
vector-con^nent  grids  is  as  follows; 


cos6 


^0*  eo»»J.+  mmh 


(la) 


sin6  =  -  cos*0 


(lb) 


(Ic) 


r* 


+  2hja^  -  (a*  -  b*)  sin* A. 


-  (a*  -b*)  lia^X 
•2-(a2-b2)na^X 


(Id) 
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<!>**<(> 


(Ic) 


y^ere,  the  parameters  a  and  b  are  the  semi-major  and  semi-minor  axes  of  the  ellipsoid  to  which 
the  geodetic  coordinates  are  referenced.  In  the  case  of  the  Juan  de  Fuca  survey,  the  reference 
ellipsoid  was  that  corresponding  to  the  World  Geodetic  System  1972  (WGS-72).  The 
corresponding  nu^etic-field  vector-component  rotation  from  geodetic  to  spherical  coordinates 
is: 


Br(r,6,(p)  =  smaB,(X,(p,h)  +  cosaBz(X,<p,h) 
Be(r,6,(p)  =  ~  cosaBzCX^cp.h)  +  sinaBz(^({>,h) 
B,(r,e,«p)  »  By(X,<p,h) 

where  the  rotation  angle  a  is  defined  as: 

a  =  X,  +6  -  Y 


(2a) 

(2b) 

(2c) 


(3) 


and  udiere  6  is  conqtuted  from  eq.  (la).  Equivalently,  the  angle  a  may  be  defined  by  the 
following  convenient  relations: 


sina  as 


~b*)  ri»X  totX 


(4a) 


¥+  cot*X+b*  *w^X 

cosa  * - r - 


(4b) 
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The  subsequent  transfonmtion  from  spherical  coordinates  (r.O.ip)  to  rectangular  coordinates 
(x,y^)  assumes  that  there  is  a  plane  tai^ent  to  the  sphere  with  its  origin  at  the  point  (R,0q,<Pq), 
^ere  R  is  the  mean  radius  of  the  Earth  as  determined  by  the  parameters  of  the  reference 
ellipsoid,  while  6o  and  <Po  define  the  center  of  the  survey  area  in  spherical  coordinates. 

The  rectangular  coordinate-system  axes  are  oriented  with  respect  to  the  spherical  l  .mate 
system  such  that  the  X-axis  points  north,  in  the  direction  of  -6,  the  Y-axis  points  east  in  the 
direction  of  +9,  and  the  Z-axis  points  vertically  down  in  the  direction  of  -r,  thus  making  a 
right-handed  system.  Consequently,  the  transformation  from  spherical  coordinates  to 
rectangular  coordinates  takes  the  following  mathematical  form: 


X  =  rsinV  cost 

(5a) 

y  s=  rsinH*  sinT 

(5b) 

z  =  R  -  rcosV 

(5c) 

udiere: 

^  s  cos~'[cos6o  cos6  +  sindo  sin6  cos((p 

^  L  «i'  J 

-  <Po)] 

(6a) 

(6b) 

and  where; 

cosT  ss  +  ^1  -  sin^T 

Irl  s  90* 

(7a) 
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cotT 


iTl  >  +  90® 


(7b) 


■  -  Jl  -  sm*T 

The  correqKmdiiig  rotation  of  the  inagnetic<fieid  vector  con^Kments  is: 

(8a) 
(8b) 
(8c) 


Bi(x,y,z)  =  -  B,  cos6o  cos<po  -  B,  cosQo  sinipo  +  B,  sinOo 


By(x,y,z)  s  -  B.  sin<po  +  ByCOsq>o 


/  I  • 

B»(x,y,z)  =  -  B,  sinBo  cos<(>o  -  By  sinBo  sin(po  -  B*  cosBo 


where: 

B,  »  B»  cosB  C0S9  -  B,  sin<p  +  B,  sinB  cos«p  (9a) 

By  s  B«  cosB  sin<|>  +  B^  cos9  +  B  sinB  sintp  (9b) 

Bx  =  -  B#  sinB  +  B,  cosB  (9c) 

After  perfoiming  these  transformations,  the  resulting  gridded  magnetic-field  components 
are  suitable  for  rectangular-harmoiuc  analysis. 


3.  THE  POTENTIAL  FIELD  MODEL 

The  three  residual  grids  corresponding  to  the  three  magnetic-field  components  of  the  Juan  de 
Fttca  area  are  presumed  to  be  related  to  each  other  through  a  single  magnetic  potential  V(x,y,z). 
This  potential  is  con^tosed  of  two  parts,  a  r^onal  potential  V|^(x,y,z)  and  a  local  potential 
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VL(x,y^)*othat: 


V(x,y,z)  «  VE(x,y,z)  +  VL(x,y,z)  (10) 

The  regional  potential  V|i(x,y,z)  includes  magnetic-field  contributions  generated  in  the  Earth's 
fluid  core  and  contributions  generated  by  induced  and  remnant  m^etization  originating  in  the 
deqper  crust  above  the  Curie  depth.  The  regional  magnetic  potential  has  characteristic 
wavelengths  that  are  longer  than  the  dimensions  of  the  survey  area.  The  local  potential  Vj^(x,y^) 
on  the  other  hand  is  presumed  to  be  due  to  Main  field  induced  and  remnant  magnetizations 
originating  in  the  Earth's  upper  crust  and  has  characteristic  wavelengths  on  the  order  of  or 
shorto'  than  the  dimensions  of  the  survey  area,  which  is  taken  to  be  rectangular  with  dimensions 
L,  x  Ly  .  Fran  these  lengths,  a  single  characteristic  dimension  L  of  the  survey  area  can  be 
defined  as: 

L «  jul;  (11) 

In  the  case  of  the  Juan  de  Fuca  survey,  L,  =  664.466  kilometers,  while  Ly  =  428.406  kilometers 
so  that  L  »  533.537  kilometers. 

3.1  The  Regional  Magnetic-Field  Model 

The  regional  magnetic-field  potential  is  taken  to  be  of  the  following  quadratic  form: 

Va(x,y,z)  »  ajix^  +  -^bnyx^x"  p,v  =  1,2,3  •  (12) 
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\«^iere  tensor  natation  has  been  used  in  v^ch  the  co<»dinates  (x,y^)  coire^nd  to  (x*  and 
\<diere  Einstein  summatioi  notadon  is  used  so  that  repeated  indices  are  to  be  summed  over. 
Thus,  for  instance: 

anx"  •  ai x'  +  ai X*  +  ai X*  =  ai x  +  aiy  +  aj z  (13) 

Using  the  single  «lgi»h  (/)  notation  (Adler,  Bazin,  and  SchifTer  [1975])  to  denote  partial 
differentiation,  so  that  for  instance  (  /p  ■  ^  ),  the  regional  magnetic  field  B,i(x,y,z),  which  is 
the  negative  gradient  of  the  regional  potential,  may  be  computed  as: 

Bi»(x,y,2)  *  -VK/,i(x,y,2)  =  -a„  -  b»»x''  (14) 

It  thus  becomes  clear  that  the  constants  -a^  are  the  regional  mean  values  of  the  corresponding 
magnetic-fidd  components.  The  local  magnetic  fidd  will  dierefore  have  a  zero  mean  for  each 
of  its  magnetic  conqionents.  Subsequently,  taking  the  positive  gradient  of  the  regiond 
nu^nedc  field  yields  the  regional  magnetic  gradient  tensor. 

^*iiv(x,y,z)  *  Bi,i/v{x,y,z)  =  -b„»  (15) 

Thus,  it  is  clear  that  the  nine  constant  coeffidents-b^,  correspond  to  the  average  gradient  or 
slope  over  the  entire  survey  area  of  each  of  the  three  magnetic  conq)onents  in  each  of  the  three 
cotffdinate  directions. 

Since  the  regional  magnetic  gradient  is  obtained  by  taking  two  gradioits  of  the  regiond 
magnetic  potentid  aid  since  the  order  in  udiich  these  gradients  are  performed  is  of  no 
conseqnence,  the  regional  magnetic>gradiait  tensor  must  be  ^mmetric.  Also,  since  the  regiond 
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magnetic  potential  miut  satisfy  Li^lace's  equadtm,  the  regional  magnetic-gradient  tensw  must 
also  have  zero  trace.  So,  this  tensor  has  only  five  independent  components  and  the  coefficients 
satisfy  the  following  syimnetry  conditions: 


b|iv  *  bv|i 

(16a) 

b^"  =  0 

(I6b) 

Experimentally,  due  to  noise  and  other  expmmental  factors,  these  relations  are  neva-  exactly 
satisfied.  The  values  obtained  via  least-squares  fitting  for  the  regional  magnetic  biases  a,^  and 
the  regional  m^etic  gradients  b^y  for  the  Juan  de  Fuca  area  are  listed  in  Table  2.  Note  that  the 
vertical  gradients  were  not  measured  since  this  would  have  required  a  second  identical  survey  at 
a  second  altitude.  However,  these  gradients  may  be  inferred  from  the  information  in  Table  2 
usiiig  die  symmetry  conditions  above.  Note  also  that  due  to  noise  and  other  measurement  errors, 
the  b,  2  coeffidem  is  not  identical  to  bj ,,  so  that  eq.  (16a)  is  not  completely  satisfied.  However, 
the  average  of  these  two  quantities  may  be  taken  as  the  best  estimate  for  both.  Table  3  is  the 
sqqinqniate  revision  of  Table  2,  which  does  satisfy  all  of  the  symmetry  relations  for  the 
regional  fidd. 

The  regional  magnetic  fidd,  as  determined  from  eq.  (12)  and  the  coeffidents  in 
Table  2,  was  removed  from  the  three  gridded  magnetic-field  component  data  sets  that 
magnetically  characterize  the  Juan  de  Fuca/Explorer  plate  regions.  The  resulting  residual  grids 
dien  define  the  local  magnetic  bdiavior  of  the  combined  regions. 
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Table  2.  Refioaal  Macoctic  Model  for  the  Nortbov  Juan  de  Fuca/Explorer  Plates  H 

Panoxieter 

Units 

Value  • 

L, 

km 

664.4660  1 

L, 

km 

428.4060  J 

L 

km 

533.5370  ^ 

z. 

km 

5.3890  1 

nT 

-32.3719  J 

^2 

nT 

14.2246  n 

nT 

14.1783  1 

bn 

nT/km 

0.1401  J 

t>0 

nT/km 

0.1589  J 

b,, 

nT/km 

—  1 

b2« 

nT/km 

-0.0280  J 

bjj 

nT/km 

-0.0509  J 

ba 

nT/km 

—  1 

b„ 

nT/km 

-0.0261  1 

bn 

nT/km 

-0.1225  J 

b„ 

nT/km 

-  1 
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Table  3.  Regional  Magnetic  Model  for  the  Northern  Juan  de  Fuca/Ezplorer  Plates 
Revision  of  Table  2  Based  on  Symmetry  Consktoratkins 


Parameter 

Units 

Value 

K 

km 

664.4660 

h 

km 

428.4060 

L 

km 

533.5370 

km 

5.3890 

*1 

nT 

-32.3719 

h 

nT 

14.2246 

aj 

nT 

14.1783 

b„ 

nT/km 

0.1401 

b„ 

nT/km 

0.0654 

bti 

nT/km 

-0.0261 

b:. 

nT/km 

0.0654 

b« 

nT/km 

-0.0509 

b^ 

nT/km 

-0.1225 

bsi 

nT/km 

-0.0261 

b,: 

nT/km 

-0.1225 

bss 

nT/km 

-0.0892 

3.2  The  Local  Magnetic-FidldModd 


It  is  now  assumed  that  the  regional  trends  and  biases  have  been  removed  from  the  gridded 
survey  The  local  potential-field  model  is  then  taken  to  be  the  solution  of  Laplace's 
in  rectangular  coordinates  and  may  be  written  in  the  following  form: 

N  M 

VL(x.y,2)  «  ss  [A*  cosn0(x)  cosmqKy)  +  B»i  co5n0(x)  sinmepiy)  + 

a>l  «-l 

(17) 

D«i  sinn0(x)  cosmqKy)  +  E»smne(x)  sinnwiKylle^-* 
where,  for  a  rectangular  area  of  dimensions  L,  x  : 


0(X)  =  ^ 

X  ^  +-^ 

(18a) 

<p(y)«^ 

(18b) 

and  whoe  the  wavenumber  k,,  is  given  by  the  relation: 

k-  *  2k  I^nSN;  ISm^M  (19) 

The  parameters  N  and  M  were  both  assigned  the  value  of  128  in  anticipation  of  the  subsequent 
use  of  die  2-dimen8ional  FFT,  while  the  parameter  L  was  previously  defined  by  eq.  (1 1). 

The  olqective  is  to  determine  the  rectangular-harmonic  coefficients  ,  D,,  , 

and  E, .  ,  which  are  real-valued  constants  with  units  of  nanoTeslas  (nT),  using  the  ^dded, 
lendual  magnedc-fidd  vector  conqKments  doived  from  the  survey  data.  These  m^etic-fidd 
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cou^QMiits  ate  related  to  the  negative  gradients  of  the  local  magnetic  potential  defined  by  eq. 
(17),  so  that: 


BL(x,y,z)  =  -VVL(x,y,2)  (20) 


Consequently,  the  three  residual  magnetic*field  conq>onents  may  be  expressed  as  follows: 


N  M 

BL.(x,y,2)  «  +  ZS  n[AMi  sinnd(x)  cosmqKy)  +  Bn  smnd(x)  sinnupiy)  - 

Bs|  B>l 

(21a) 

Dob  cosn6(x)  cosm(p(y)  -  En  cosn6(x)  sinm(p(y)]e'‘*”‘ 


N  M 

BLy(x,y,2)  «  +  EX  ni[A«a  cosn0(x)  sinmiiKy)  -  Bn  co$d6(x)  cosnup(y)  -i- 

■a|  Ba| 

(21b) 

Dn  sinn0(x)  sinmqKy)  -  En  sinn0(x)  cosm(p(y)]e^“* 


N  N 

Bi,(x,y,z)  =-•^22  1^[Abd  cosn6(x)  cosmqKy)  +  Bn  cosnOCx)  sinm(p(y)  + 

Ba|  Bv| 

(21c) 

Dn  sinne(x)  cosm(p(y)  +  En  sinn0(x)  sinm((Ky)]e'‘-* 


For  convoiience  in  the  following  discussion,  the  subscript  L  will  be  suppressed,  and  it  will 
be  assumed  from  this  point  onward  that  we  are  concerned  with  the  local  magnetic  potential  and 
the  local  fields  and  gradients  that  can  be  derived  from  it.  We  then  proceed  by  noting  that  the 
rectai^ular-haimcMiic  coefficients  of  the  local  magnetic  potential  can  be  conq>uted  from  any  one 
of  the  three  gridded  conqxments  of  the  geomagnetic  residual  field,  which  are  presumed  to  have 
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been  niMasuied  at  stnne  neariy  constant  survey  altitude,  H,  (i.e.,  z  =  -H)  above  the  x-y  plane.  If 
the  X*conq)onent  of  the  residual  magnetic  field  is  used,  then,  using  the  usual  orthogonality 
conditions  for  the  sine  and  cosine  functions,  these  coefficients  are  determined  as  follows: 


J  /  cosm(pB,(0,<p,-H)ded(p 

(22a) 

J  J  sinn6  sinnup  B.(6,q>,-H)d0d(p 

(22b) 

“  ”  n^T  ***"”  J  J  cosn6  cosoKp  B*(0,<p,-H)d0d(p 

(22c) 

~  “  7^  *"'*"**  J  /  cosn0  sinm{pBi,(0,<p,-H)d0d(p 

(22d) 

Alternatively,  using  the  Y>component  magnetic  residuals,  these  coefficients  are  determined  as: 

L  .  f*  f* 

ro  Jx  J  J  sinro9By(0,(|>,“H)d0d(p 

(23a) 

L  f*  f* 

"  “  rai^  J  J  cosm(|)By(0,(p,-H)d0d(p 

(23b) 

J*  J  sinn0  smm(pBy(0,(|>,-H)d0d(p 

(23c) 

E—  =  -  — J  J  suin0  cosm9By(0,(p,-H)d0d(|> 

(23d) 

The  diird  alternative  uses  the  Z-component  of  the  residual  magnetic  field  and  yields  the 
following  egqnDessions  for  these  coefficients: 
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A»i  =  - 


B«.  =  - 


^gk«iH  J  J  cosn0  cosm<pBt(0,<p,-H)d0d(() 
J  J  cosn0  suiffl(pBi(0,(p,-H)d0d(|> 


D_  - 2 — 

“  llLkn 


sinn0  cosin(pBz(0,q>,-H)d0d(p 


lolzr  I  J  sinn0  smm<pB2(0,(p,-H)d0d9 


(24a) 


(24b) 


(24c) 


(24d) 


Generally,  the  magnetic  residual  component  exhibiting  the  largest  an^litude  is  the  preferred  one 
for  computing  the  rectangular-harmonic  coefficients.  This  choice  will  yield  smaller  coefficient 
errors  due  to  the  smaller  percentage  of  noise  intrinsic  to  the  data  for  that  magnetic  component 
relative  to  the  others.  The  Z-component  grid  was  chosen  for  this  purpose  in  the  Juan  de 
Fuca/Explorer  region.  The  resulting  coefficients  are  partially  listed  to  degree  15  in  Table  4. 

Evaluation  of  the  above  expressions  for  the  rectangular-harmonic  coefficioits  is 
accomplished  with  the  aid  of  the  2-dimensional  FFT.  First,  note  that  each  residual, 
magnetic-coiiq>onent  grid  can  be  represented  in  the  form  of  a  2-dimaisional,  ffiiite  Fourier  series 
as  follows: 


N  M 


Bz(0,q>,-H)  =  Z  2  bxv»(-H)e‘<''»*'‘^ 


v»-N  |ts-M 


N  M 

B,(0.<p,-H)  =  2  Z  b,vH(-H)e‘<''**'‘*» 

v«-N  |i>-M 


(25a) 


(25b) 
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Table  4.  Rectangiilar>Harmoiiic  Model  Coefflcients  of  the  Juan  de  Fuca/Explorer  Plate 
Regions,  Low  Degree  and  Order  Terms  Only;  Units:  nanoTeslas 


n 

m 

A.. 

B.. 

D- 

E.. 

1 

1 

-6.50182 

0.90794 

10.10035 

6.98122 

1 

2 

-4.20178 

-4.19915 

-8.77588 

1.38671 

1 

3 

-0.17999 

2.12747 

2.20755 

2.62100 

1 

4 

-0.26331 

-3.23829 

0.96876 

-1.44732 

1 

5 

-2.97880 

-0.39970 

-1.28883 

-1.16297 

1 

6 

0.24739 

1.26363 

-1.31373 

0.99827 

1 

7 

0.07261 

-0.82978 

0.78405 

0.23625 

1 

8 

-0.19026 

0.76457 

-0.15335 

0.17747 

1 

9 

0.16159 

-0.28253 

0.38965 

-0.44789 

1 

10 

-0.51639 

0.44997 

-0.63982 

-0.07680 

1 

11 

0.16159 

0.12405 

0.04508 

-0.06916 

1 

12 

-0.51639 

0.10945 

-0.34906 

-0.34367 

1 

13 

0.05547 

0.09868 

-0.18891 

0.07211 

1 

14 

-0.32374 

0.08613 

-0.04772 

0.05725 

1 

15 

-0.03293 

-0.11430 

-0.03724 

0.07085 

2 

1 

15.71674 

17.28025 

18.39594 

-5.71595 

2 

2 

-2.36503 

-4.43721 

-0.32715 

7.12863 

2 

3 

-1.61875 

-0.47278 

-0.32715 

-4.94449 

2 

4 

-0.99136 

-1.33616 

0.32225 

2.52737 
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Table  4.  Rectangular-Harmonic  Model  Coeflicienta  (con.) 


2 

5 

-0.75027 

-0.20994 

-2.13968 

0.08772 

2 

6 

-0.53263 

-0.76496 

-0.98668 

-0.80986 

2 

7 

-0.60334 

0.69700 

-0.08925 

0.71491 

2 

8 

0.22274 

0.13464 

0.32046 

-0.32020 

2 

9 

0.06444 

-0.35671 

-0.00394 

0.16452 

2 

10 

0.00910 

-0.42301 

0.58391 

-0.04872 

2 

11 

-0.23317 

-0.44016 

0.12173 

-0.23388 

2 

12 

0.03802 

-0.16201 

0.19692 

0.22992 

2 

13 

-0.04182 

-0.24281 

0.26783 

-0.10191 

2 

14 

-0.20038 

0.05171 

-0.05880 

-0.09223 

2 

15 

-0.02496 

0.16278 

0.00662 

0.04175 

3 

1 

14.69034 

6.64214 

26.20008 

-6.60164 

3 

2 

-6.92724 

-1.22531 

0.75517 

-3.78142 

3 

3 

-1.85056 

-1.13347 

-0.55870 

0.53041 

3 

4 

-0.58444 

0.98213 

-0.29600 

-1.29544 

3 

5 

0.07912 

1.05613 

-1.77336 

-0.00432 

3 

6 

0.39846 

0.04386 

0.10321 

1.21775 

3 

7 

-0.08147 

-0.52859 

0.56867 

0.05593 

3 

8 

-0.11638 

-0.16341 

0.51915 

-0.13795 

3 

9 

-0.39917 

0.27863 

-0.31902 

-0.17071 

3 

10 

0.41564 

0.78612 

-0.43821 

0.49178 
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3 

11 

0.50470 

0.05392 

0.02390 

0.44870 

3 

12 

0.07068 

0.23216 

-0.07872 

0.03265 

3 

13 

0.20408 

0.00148 

-0.11725 

0.17782 

3 

14 

-0.05528 

-0.08068 

0.24916 

0.06604 

3 

15 

0.03076 

0.05477 

-0.00173 

-0.22982 

4 

1 

-6.72438 

5.38535 

2.12106 

-10.24770 

4 

2 

-1.76723 

2.90900 

7.01444 

0.53800 

4 

3 

-0.40940 

0.64107 

-0.80526 

-0.26738 

4 

4 

-0.65961 

-1.39914 

-0.59056 

-0.91643 

4 

5 

-0.88687 

0.32765 

-0.42612 

-0.86223 

4 

6 

0.43625 

-0.16263 

-0.17006 

-0.43458 

4 

7 

-0.66605 

0.63265 

-0.45634 

0.21541 

4 

8 

0.65868 

0.08300 

0.04932 

0.22266 

4 

9 

-0.35682 

-0.32375 

0.06753 

0.00346 

4 

10 

-0.31478 

0.00122 

0.24568 

-0.32948 

4 

11 

-0.22759 

0.11451 

-0.23232 

-0.43214 

4 

12 

0.05813 

-0.13130 

0.02407 

-0.02362 

4 

13 

-0.18008 

-0.01142 

-0.00368 

-0.20583 

4 

14 

0.15060 

0.10593 

-0.10623 

0.06501 

4 

15 

-0.02681 

-0.15914 

0.02728 

-0.00179 

5 

1 

-9.30405 

-2.29585 

0.19474 

-11.34965 
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5 

2 

3.86618 

1.33161 

0.43339 

-1.80079 

5 

3 

-1.41536 

1.66478 

-0.85717 

1.08008 

5 

4 

1.02127 

-1.91231 

0.59884 

1.66993 

5  , 

5 

-1.18006 

-0.42611 

0.38220 

-0.17926 

5 

6 

0.42373 

-0.07830 

-0.19304 

0.17104 

5 

7 

-0.71082 

-0.53879 

0.81771 

-0.81005 

5 

8 

-0.32892 

0.48116 

-0.34766 

-0.26324 

5 

9 

0.04880 

0.48425 

-0.18618 

0.12710 

5 

10 

0.66700 

0.13226 

-0.01505 

0.66169 

5 

11 

0.06542 

-0.41894 

0.38528 

0.12194 

5 

12 

-0.10911 

0.17215 

-0.08716 

0.02790 

5 

13 

0.12715 

0.01377 

0.00680 

-0.00514 

5 

14 

-0.02725 

-0.00738 

-0.12867 

0.05027 

5 

15 

0.08473 

0.01906 

0.06906 

0.12909 

6 

1 

-1.28958 

-1.31317 

-1.35462 

-0.85454 

6 

2 

2.25894 

-2.98705 

-1.59060 

2.22455 

6 

3 

-1.39635 

2.90327 

2.18386 

-2.61767 

6 

4 

1.09335 

0.01516 

-1.91373 

1.55669 

6 

5 

-0.43473 

0.33791 

0.72297 

0.51169 

6 

6 

0.02327 

-0.55756 

-0.15552 

-0.28759 

6 

7 

0.03012 

0.38200 

-0.19865 

0.46374 
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6 

8 

0.67396 

-0.27972 

0.42758 

0.19033 

6 

9 

-0.40740 

-0.34895 

0.14588 

-0.70268 

6 

10 

-0.23443 

0.20503 

-0.43304 

-0.22173 

6 

11 

-0.12858 

0.19569 

-0.12127 

-0.05647 

6 

12 

0.04700 

-0.03271 

0.02773 

0.01470 

6 

13 

-0.13390 

-0.00309 

0.036  72 

-0.17762 

6 

14 

-0.02297 

0.09231 

-0.10085 

0.00624 

6 

15 

-0.00697 

0.11875 

-0.03650 

-0.05949 

7 

1 

0.42043 

1.76266 

1.61781 

-0.70011 

7 

2 

-1.73306 

-0.44276 

-2.84195 

0.38552 

7 

3 

1.58631 

-0.29854 

0.83468 

-1.85681 

7 

4 

-0.70798 

-0.00683 

0.01928 

-0.12694 

7 

5 

-0.84105 

0.81417 

-0.10592 

-0.70527 

7 

6 

0.50023 

-0.00103 

-0.18092 

0.91178 

7 

7 

-0.70702 

-0.11124 

0.19987 

-0.58793 

7 

8 

0.16677 

0.56344 

-0.66850 

0.22586 

7 

9 

0.20803 

-0.31914 

0.35395 

0.46763 

7 

10 

0.06335 

-0.01063 

0.24076 

0.11647 

7 

11 

-0.04250 

-0.28750 

0.24752 

-0.11983 

7 

12 

-0.13971 

0.28525 

-0.19586 

0.03015 

7 

13 

0.27651 

-0.06010 

0.06477 

0.17456 
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7 

14 

-0.06703 

-0.16502 

0.14928 

-0.00955 

7 

15 

-0.03639 

-0.08032 

0.07791 

-0.17341 

8 

1 

3.31603 

5.58220 

-0.44074 

-0.79490 

8 

2 

-0.95824 

-4.46747 

1.86770 

-0.39213 

8 

3 

-1.34737 

0.91580 

-1.44776 

0.25701 

8 

4 

1.10121 

-0.29209 

0.70388 

-0.42378 

8 

5 

-0.30687 

-0.29736 

0.50434 

-0.69129 

8 

6 

-0.13969 

-0.04296 

-0.53447 

-0.26891 

8 

7 

-0.03712 

0.36213 

0.02546 

0.62256 

8 

8 

0.25247 

-0.45761 

0.44310 

-0.11620 

8 

9 

-0.34616 

0.32148 

-0.11331 

-0.26649 

8 

10 

0.18843 

0.23557 

-0.17133 

-0.20744 

8 

11 

0.10870 

0.05990 

-0.11914 

0.13291 

8 

12 

-0.02481 

-0.19952 

0.23894 

-0.17592 

8 

13 

-0.19212 

0.19060 

-0.23797 

-0.12452 

8 

14 

0.23414 

0.06747 

-0.04150 

0.19298 

8 

IS 

0.01644 

-0.09467 

-0.00029 

0.06760 

9 

1 

-2.18537 

3.48995 

-1.57473 

-4.78597 

9 

2 

2.16019 

-2.49957 

0.41614 

2.09978 

9 

3 

-2.26432 

1.03862 

0.00126 

-1.16616 

9 

4 

0.77841 

-0.74632 

-0.37452 

0.35329 
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9 

5 

-0.17226 

0.36005 

-0.29210 

-0.29683 

9 

6 

0.30127 

-0.45394 

0.18381 

0.83077 

9 

7 

-0.83058 

0.26026 

0.17586 

-0.87825 

9 

8 

0.58998 

0.40076 

-0.69287 

0.66965 

9 

9 

-0.17553 

-0.44890 

0.51060 

-0.08050 

9 

10 

0.04407 

0.08372 

-0.09538 

-0.07233 

9 

11 

0.03944 

-0.07694 

-0.04026 

-0.11101 

9 

12 

0.19002 

0.01810 

-0.12889 

0.32356 

9 

13 

-0.04717 

-0.21440 

0.33878 

-0.09185 

9 

14 

-0.16804 

0.06191 

-0.16683 

-0.13568 

9 

15 

-0.03636 

0.08599 

-0.00703 

0.00700 

10 

1 

-1.53264 

0.88371 

1.29867 

-0.03071 

10 

2 

1.50872 

-0.92994 

0.03802 

-0.20592 

10 

3 

-1.81144 

-0.57259 

0.06813 

-0.74301 

10 

4 

0.16088 

0.55026 

-0.04285 

-0.86953 

10 

5 

-0J0456 

-0.23105 

0.07705 

-0.90593 

10 

6 

0.37293 

0.55645 

-0.49604 

0.22612 

10 

7 

0.07846 

-0.58866 

0.17024 

0.33061 

10 

8 

-0.52099 

0.28158 

-0.04822 

-0.55390 

10 

9 

0.14498 

0.07277 

-0.16160 

0.15475 

10 

10 

-0.02714 

0.06984 

-0.10329 

-0.01712 
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10 

11 

•0.04698 

-0.10054 

-0.00504 

0.13541 

10 

12 

-0.18429 

0.17042 

-0.07670 

-0.12646 

10 

13 

0.25393 

0.06462 

-0.05330 

0.12606 

10 

14 

-0.02404 

-0.12619 

0.11820 

0.00188 

10 

15 

-0.00626 

0.00404 

-0.00203 

-0.00244 

11 

1 

0.68794 

2.00314 

-1.82657 

1.29417 

11 

2 

0.68014 

-2.00106 

0.49140 

-0.01759 

11 

3 

-0.64070 

-0.22335 

0.41249 

-0.07980 

11 

4 

0.40475 

0.12301 

-0.11870 

-0.20360 

11 

5 

-0.43983 

-0.53327 

-0.08233 

-0.20360 

11 

6 

0.04881 

0.10360 

0.04355 

0.38457 

11 

7 

0.13693 

0.20792 

-0.26172 

-0.22968 

11 

8 

0.19715 

-0.11042 

0.40007 

0.01963 

11 

9 

-0.21894 

-0.14376 

0.01066 

0.17168 

11 

10 

-0.08586 

0.18476 

-0.05054 

-0.21125 

11 

11 

0.22761 

-0.10226 

0.16959 

0.12777 

11 

12 

-0.12130 

-0.10635 

0.08399 

0.06518 

11 

13 

0.01405 

0.12173 

-0.14097 

-0.09379 

11 

14 

0.04797 

-0.04983 

-0.02296 

0.06785 

11 

IS 

-0.02768 

0.07813 

-0.00798 

-0.01461 

12 

1 

-1.90696 

2.52971 

-4.09107 

3.08946 
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12 

2 

1.59213 

-0.73171 

1.46118 

0.80026 

12 

3 

-0.14431 

-0.39781 

-0.69909 

0.13509 

12 

4 

-0.21148 

-0.33016 

-0.16375 

0.21105 

12 

5 

-0.36754 

-0.57927 

0.13807 

0.37689 

12 

6 

0.33746 

0.01655 

0.01266 

-0.07096 

12 

7 

-0.19403 

-0.23841 

-0.15226 

-0.02569 

12 

8 

-0.20600 

0.24767 

0.03850 

-0.03410 

12 

9 

0.14263 

-0.13189 

0.02403 

0.04626 

12 

10 

-0.15607 

0.04710 

0.10123 

-0.07676 

12 

11 

0.10118 

0.10544 

-0.21553 

0.00513 

12 

12 

-0.01595 

-0.14016 

0.04774 

0.05590 

12 

13 

-0.08750 

0.05428 

-0.05250 

-0.06875 

12 

14 

-0.07930 

-0.01179 

-0.01840 

0.07531 

12 

15 

-0.05007 

-0.02926 

0.04668 

-0.07477 

13 

1 

-4.67324 

4.92206 

-0.48776 

1.65893 

13 

2 

2.04069 

-0.47986 

-0.02500 

-0.68042 

13 

3 

-1.38261 

0.21989 

-0.59911 

0.40972 

13 

4 

0.21093 

0.02153 

0.64767 

0.04682 

13 

5 

-0.49765 

0.15527 

-0.26342 

-0.33375 

13 

6 

0.11866 

0.28142 

0.20178 

0.04789 

13 

7 

-0.12871 

-0.11888 

-0.28160 

-0.11622 

30 


Table  4.  Rectangnlar-Haraoik  Model  Co^nctents  (con.) 


13 

8 

0.07702 

0.09139 

-0.03837 

0.21612 

13 

9 

-0.02636 

-0.06746 

-0.03337 

-0.16914 

13 

10 

0.10606 

0.11825 

0.03435 

0.12884 

13 

11 

-0.01994 

-0.11033 

0.06455 

-0.12492 

13 

12 

-0.08246 

0.07670 

-0.05334 

0.12445 

13 

13 

0.05388 

0.01529 

0.06806 

-0.03654 

13 

14 

-0.03364 

-0.03445 

-0.08154 

0.04046 

13 

15 

-0.00051 

-0.02318 

0.07854 

0.00493 

14 

1 

-2.1U5y 

5.23947 

2.84025 

-1.75222 

14 

2 

2.05113 

-1.83951 

-1.14319 

-0.07571 

14 

3 

-1.51667 

-0.12111 

1.02495 

0.36313 

14 

4 

0,75293 

0.22577 

0.05374 

-0.51673 

14 

5 

-0.49578 

-0.32312 

-0.04813 

0.14631 

14 

6 

0.37889 

0.15998 

0.29040 

-0.16016 

14 

7 

-0.34865 

-0.31172 

-0.06350 

-0.06328 

14 

8 

-0.05944 

0.34002 

-0.07941 

0.01225 

14 

9 

0.10337 

-0.23342 

0.11621 

0.03923 

14 

10 

-0.15907 

0.09641 

-0.12157 

-0.24132 

14 

11 

0.09470 

0.00729 

-0.07918 

0.15628 

14 

12 

-0.06836 

0.00526 

0.13284 

-0.10855 

14  13 


0.02S88 


0.04462 


-0.14339 


-0.02591 
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14 

14 

0.00271 

0.00534 

0.06132 

0.10086 

14 

15 

0.01955 

-0.01534 

0.02426 

-0.01894 

15 

1 

-0.49015 

1.55673 

2.55386 

-5.02614 

15 

2 

0.11961 

-0.83393 

-2.46010 

0.95284 

15 

3 

-0.15853 

0.47131 

1.06937 

0.71696 

15 

4 

0.38033 

-0.33976 

-0.27622 

-0.51193 

15 

5 

-0.51728 

-0.04921 

-0.04423 

0.35401 

15 

6 

0.34257 

0.22986 

-0.06258 

0.05729 

15 

7 

-0.01659 

-0.26012 

0.16216 

0.07511 

15 

8 

-0.17459 

0.09203 

0.07409 

-0.04397 

15 

9 

0.11633 

0.02952 

-0.03338 

0.05132 

15 

10 

-0.06943 

-0.10939 

0.08129 

0.00419 

15 

11 

0.05010 

-0.10604 

0.00194 

-0.01563 

15 

12 

0.06284 

-0.10455 

0.05631 

0.01565 

15 

13 

-0.11932 

0.00911 

0.00218 

-0.07397 

15 

14 

0.04238 

0.04580 

-0.04159 

0.00744 

15 

15 

-0.01570 

-0.01475 

0.00892 

-0.00096 
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V^N  ^  m~U 


(25c) 


The  complex-valued  Fourier  coefficients  b,  ^  ^(-H),  by  ^  y(-H),  and  b^  ha’*^c 

units  of  nanoTeslas  (nT),  are  known  quantities  since  they  can  be  obtained  by  taking  the  FFT  of 
their  req)ective  X-,  Y-  or  Z-c(mq)onent,  magnetic-residual,  survey-generated  grids. 
C(»isequently,  \riien  eqs.  (25a,  25b,  and  25c)  are  inserted  into  eqs.  (22a  through  22d,  23a 
through  23d,  and  24a  through  24d),  the  rectangular-harmonic  coefficients  can  be  analytically 
evaluated  in  terms  of  their  re^)ective  Fourier  coefficients  by  converting  the  sine  and  cosine 
functions  in  the  above  expressions  to  their  exponential  forms  and  enq>loying  the  Dirac  delta 
function  identity: 


J*  =  2«5(X-X') 


(26) 


over  all  integrals  involving  the  parameters  6  and  9.  The  integrals  to  be  evaluated  are  of  just  the 
following  two  basic  forms: 


J  e~'*^  sinicAdQ  =  in6(K  +  X)  -  iJtS(K  -  X) 
J  cosxDdn  =  k6(k  +  X)  +  x6(k -  X) 


Ctmsequently,  eqs.  (22a)  through  (22d)  are  evaluated  as: 


(27a) 


(27b) 
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Am  *  + 


"S’  ~ 


+  bii,^a  biAB  6^ 


(28a) 


“  +  "^^bx,-«.-«  —  b|,s,^  —  bii.-«ai  +  b«,««^e''*"® 


(28b) 


)m  =  +  +  K».-m  +  b,.-,..  +  b,.,,  je'^ 


(28c) 


i-  =  -  ■x(b..-^^  +  b,.,.-.  -  b,.-...  -  b....]( 


(28d) 


Equations  (23a)  through  (23d)  become: 


(29a) 


+  by.*.-«  +  by.^,  +  by.^»je'‘-“ 


(29b) 


(29c) 


•“  by.«,-«  +  by,_ 


-  by..,»je'‘-“ 


(29d) 


The  evaluation  of  eqs.  (24a)  through  (24d)  yields  the  third  alternative  means  of  conq)uting  the 
rectangular-harmonic,  nu^etic-potentiai  coefficients: 


Am  ®  ~  kML  bi.n,-Bi  +  bz,-a,«  +  C*'"** 


(30a) 


Bm  »  -i- 


S: 


+  bz.a.-a  -  bi.-^. 


(30b) 
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D-  -  +  ^  -  b^...)e'^«  (30c) 

E-  =  +  ^ -  bz.,.^  -  bz.^.  +  bz...-)e^-«  (30d) 

The  rectangular-hannonic  coefficients  of  the  local  magnetic  potential  are  thus  determined 
through  a  two-step  process; 

(1)  Perform  an  FFT  on  one  of  the  three  available  gridded  vector  components  of  the 
local  magnetic  field. 

(2)  Use  the  Fourier  coefficients  obtained  in  step  1  to  evaluate  the  rectangular-hannonic 
coefficients  of  the  local  magnetic  potential  using  the  conesponding  set  of 
equations  derived  above  for  these  coefficients. 

3.2.1  Computing  the  Local  Magnetic-Field  Vector  from  the  Model 

Having  evaluated  the  coefficients  of  the  local  rectangular-harmonic  model,  it  is  now  desirable 
to  reconq)ute  all  three  vector-magnetic  components  ,  and  B, .  The  reason  for  this  is  that 
the  grids  derived  from  the  survey  data  are  not  entirely  consistent  with  each  other  because  of 
instrument  noise  and  other  experimental  factors.  However,  the  magnetic-field  components 
derived  from  the  model  are  necessarily  consistent  with  each  other  by  definition.  As  a 
consequence,  geophysical  interpretations  of  a  region  based  on  least-squares  or  other 
inverse-modeling  techniques,  using  magnetic-field  parameters  derived  in  this  way,  are  less 
ambiguous  due  to  the  reduced  noise  level  of  the  inverted  data.  Secondly,  it  is  sometimes 
desirable  to  iq>ward  or  downward  continue  the  original  survey  data  to  some  other  altitude  H*.  A 


35 


F — 


local  potential-field  model  pennits  this  procedure  to  be  performed  with  relative  ease  using  FFT 
procedures  and  specially  designed  filters,  if  they  are  needed. 

We  now  assume  that  the  rectangular-harmonic  coefficients.  A,. ,  ,  and  E,. 

are  known  andsedc  to  determine  their  corresponding  Fourier  coefficients  by^^C-IT), 

and  bj  ^  Bcpiations  (25a,  25b,  and  25c)  can  easily  be  inverted  via  the  usual  techniques, 

which  involve  the  use  of  eq.  (26)  to  yield: 

b,vn(-H0  =  j‘  j‘  e-‘<''»^'‘*>B.(e.((>.-H')ded<p  (31a) 

b,v*(-H')  =  (4n2)-«  f*  f"  e-‘(''**'‘*>By(0,((>,-HOd0d(|>  (31b) 

J-HC  </-4| 

b,v,i(-H0  =  (4n*)-'  f*  f*  e-‘(''**'‘«B,(0.<f>,-HOd0d{p  (31c) 

J-H  J-H 

where  we  tenqKnrarily  drop  the  tensor  notation  and  explicitly  reference  parametors  related  to  the 

magnetic-vector  conqxments  by  the  subscripts  x,  y,  and  z,  while  parameters  related  to  the 

magnetic-gradient  components  will  be  explicitly  refoenced  by  the  subscripts  xx,  xy,  xz,  yz,  etc. 

Gredc  letters  will  then  be  used  to  reference  either  Fourier  coefficients  or  rectangular-harmonic 

coefficients  with  the  understanding  that  rectangular-harmonic  coefficients  with  indices  equal  to 
or  less  than  zero  are  defined  to  be  zero.  Greek  indices  now  have  the  ranges:  -N  ^  v  ^  +N 

and  -M  s;  p.  <:  +M ,  where  for  the  Juan  de  Fuca/Explorer  plate  region  N  -  M  128. 

After  inserting  eqs.  (21a,  21b,  and  21c)  into  eqs.  (31a,  31b,  and  31c),  respectively,  and 

porforming  the  necessary  integrations,  vdiich  are  of  the  following  three  basic  forms: 
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I 


J  J  coslcnd^dn  =  {50  -v)6(k-^)  +  5£(-j)-v]5(k-n) 


+  8(j  -  v)5[(-k)  -  h]  +  5CH)  -  v]5[(~k)  -  4]} 


(32a) 


J  J  cosj§  sinkn  d^di]  =  -  {5(j  -  v)5(k  -  4)  +  6[(-j)  -  v]6(k  - 


-50  -  v)8[(-k)  -  ^i]  -  6CH)  -  v]5[(-k)  -  ^i]} 


(32b) 


f  f  e-i(v5+iin)  sinkn  d^<tn  =  -  {50  -  v)5(k  -  4)  -  6C(-j)  -  v]5(k  - 

J  -« 


H) 


-50  -  v)6[(-k)  -  fi]  -  6CH)  -  v]5[(-k)  -  4]} 


(32c) 


where  again  the  Dirac  Delta  function  identity  of  eq.  (26)  has  been  used,  we  obtain  with  mininMl 
effort  the  following  results: 


biiv|i(”H  )  *  +  Av,-(i  +  +  ^Bv,|i  +  —  Bv,-|i  “ 

I, 

-  ^Dv.ii  +  D_v,|i  +  Dv.-|i  +  +  i^Ev.n  -  E-v.ii  -  Ev.-<i  +  E-v.-n  jje’*”*® 


(33a) 


bjfV|l(— H  )  —  +  ^I^Av.ll  +  A-V,|l  +  Ay,^  +  A— ~  ^Bv,|l  +  B^,|l  “  By,.^  “  B_y,^^ 

I 

+  ^Dy,,l  -  Dwy.ii  +  Dy.^  -  D-y.^  j  +  i  ^Ey.,1  -  E-y.n  -  Ey.-ji  +  E_y,_^  j  j  C"*’"**'’ 


(33b) 


bzv|l(— H^)  »  +  ■  g— '  ^Ay,||  +  A-y,||  +  Ay,^  +  A-v,'^^  —  1  ^By,||  +  B-v,|i  -■  By,-|i  —  B^,-|i^ 

—  i^Dy.ji  —  D^,||  +  Dy,^  —  —  ^Ey,|i  —  E-y,||  —  Ey,-)t  +  E-y.^]]e-^“' 


(33c) 
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Mdiere  it  is  understood  that  the  rectangular-hannonic  coefficients  are  defined  to  be  zero  if  any 
index  is  less  than  or  equal  to  zero.  The  above  result  was  obtained  using  the  fact  that  1^^^  =  ^  = 

k.  ^ 

Having  thus  determined  the  Fourier  coefficients  from  the  model,  an  iqiplication  of  the  Inverse 
FFT,  which  is  essentially  an  application  of  eqs.  (2Sa,  25b,  and  2Sc)  with  H  replaced  by  H*  (the 
new  upward/downward  continued  altitude),  yields  the  desired  self-consistoit  set  of  uniformly 
spaced  grids  for  the  local  magnetic-freld  vector  conqKments.  For  the  Juan  de  Fuca/Explorer 
region,  th^  are  illustrated  in  Figures  3a  ,  3b  ,  and  3c  ,  corresponding  to  the  X-,  Y-,  and 
Z-components  at  the  original  survey  altitude  H*  =  H. 

3.2,2  Computing  the  Magnetic-Gradient  Tensor  from  the  Model 

The  magnetic-gradient  tensor  ^x,y,z)  is  the  gradient  of  the  magnetic  field  B(x,y,z).  That  is: 

Jf(x,y,z)  =  VB(x,y,z)  (34) 

In  tensor  notation,  this  equation  becomes; 

^ii*(x,y,z)  =  B,i/»(x,y,2)  =  -  V/M/»(x,y,z)  (35) 

Since  interchanging  the  order  of  differentiation  of  the  magnetic  potential  does  not  affect  the 
result,  die  mi^netic-gtadient  tensor  is  necessarily  symmetric.  The  local  magnetic  potential  also 
must  satisfy  Ljqilace's  equation.  This  is  equivalent  to  stating  that  the  trace  of  the 
magnetic-gradient  tensor  must  be  zero.  Contracting  the  indices  p.  and  v  yields  the  trace: 
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Figare3a.  Juan  de  Fuca  X-Con^xment  Magnetic  Fidd;  Units:  nT 


Juan  de  Fuca  Y-Oxnponait  Magnetic  Field;  Units:  nT 


Component  MagiKtic  Field;  Units:  nT 


ji  *  0 


(36) 


The  symmetry  and  trace  conditions  <»  the  local  magnetic  gradient  lf(x,y^)  reduce  dw  nutid)er 
of  its  independent  con^Kments  to  five  from  the  original  nine,  just  as  was  the  case  widi  die 
tegitmal  magnetic  gradient.  Of  the  five  indqiendent  com|>onents,  two  must  be  diagonal 
elements,  and  three  must  be  off-diagonal  dements.  We  will  neverthdess  examine  all  duee 
diagonal  components  as  well  as  diree  off-diagonal  cooqionents  for  a  total  of  six  conqxments. 
Then,  taking  the  indicated  derivatives  of  the  local  magnetic  potential,  the  five  independent 
conqionents  plus  one  extra  diagonal  conqxment  take  the  following  forms; 


«»«(x,y,z) 


2  X  cosm«p(y)  +  B..  cosn6(x)  sinm^fy)  + 

(37a) 

Dm  sinn6(x)  cosmipfy)  ■¥  £,■  smn6(x)  suim(p(y)]e^* 


N  M 

^m7(x,y,z)  *  +  "^X  X  sinii6(x)  sinmtpfy)  -  Bm  sinn0(x)  cosmipCy)  - 

a«l 

(37b) 

Dm  cosn6(x)  sinmtp^)  +  Em  cosd6(x)  cosmqKy)]^^”' 


^«i(x,y,z)  =  -  X  X  “k,«[AM  sinnefx)  cosmtpfy)  +  B,.  sinn6(x)  sinmtpfy)  - 

mml 

(37c) 

D.a  oosn9(x)  cosmipfy)  -  £..  cosnOfx)  sinnMp(y)]e^-' 


dFyy(x,y,z)  *  -  -2^  J)  ni*[A,a  cosn6(x)  cosmqKy)  +  Bm  cosii6(x)  sinm<p(y)  + 

..I  ..I 

(37d) 

Daa  sinn6(x)  cosm<p(y)  +  £..  sinii6(x)  sinm(p(y)]e'^~‘ 
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^»i(x,y,2)  »■-  2  ®k««CA««  cosn0(x)  sinm(p(y)  -  B..  cosii6(x)  cosin(p(y)  + 


»■!  «■! 


Daa  sinn6(x)  smm(p(y)  -  £.■  smn6(x)  cosm<p(y)]e^*" 


(37c) 


^t.(x,y.z)  «  +  ^  kSaiCAaa  cosi^x)  cosm<p(y)  -t-  Baa  co8m6(x)  sinm(p(y)  + 


Daa  sinn6(x)  cosnupCy)  +  Eaa  sinn6(x)  sinm(p(y)]e'‘‘*' 


On  the  other  hand,  each  of  the  magnetic-gradient  components  may  be  written  in  the  form  of  a 
2-dimensi(mal  finite  Fourier  Transform  at  an  altitude  z  =  -H*  which  may  possibly  be 
tq>ward/downward  continued  from  the  survey  altitude  z,  =  -H ,  so  that: 


N  M 

«aa(e,<P,-H0  =  2  £ 

v»-N  iia-M 


N  M 
v»-N  fiv-M 


N  M 
v*-N  |i«-M 


N  M 

«,y(0,q),-HO  =  2  E  Pyyvi.(-H0e“''**'‘’> 

V«-N  lia-M 


t)  Kt 

«„(0,(|>,-H')  *  2  E  Py.v.i(-H')e‘<''**'‘»> 

v«-N  n—M 


(38a) 


(38b) 


(38c) 


(38d) 


(38e) 


N  M 

«„(e,<p,-HO  =  2  Z  (38f) 

v»-N  |1«>M 


Uang  eq.  (26),  eqs.  (37a  through  37f)  can  be  inverted  to  yidd  the  Fourier  coeffidents  in 
terms  of  die  rectangular>harmonic  coeffidents  A,. ,  B.. ,  D,. ,  and  yidding; 


P.,v,.(-H0  *  (4k^)-'  J"  e-‘<''**'‘*»48„(e,9,-H0dedip 

(39a) 

P.,v,i(-H0  *  (4tt*)-> 

(39b) 

Px.v,»(-H0  »  (471*)-'  f  r  e-‘<''**'‘««„(e,(j>.-H0dedq> 

J-n  Jhc 

(39c) 

Pyyv,i(-H0  *  (4rr*)-'  ^y,(0.9.-HOded«|) 

(39d) 

Py,y,»(-H')  s  (47t*)-'  f  r  e-‘<''**'‘»)«y*(e,(p,-H')de(k(> 

(39e) 

P«y,»(-H')  =  (471*)-'  j*  J‘  e-'<'''‘^'‘'««„(e,(|>,-H')de(k(> 

(390 

Biserting  eqs.  (36a  through  360  into  eqs.  (38a  through  380,  respectively,  and  using  eqs.  (27a 
and  27b)  to  evaluate  the  resulting  integrals,  die  Fourier  coeffidents  for  the  magnetic  gradioit 
comqponents  reduce  in  a  straightforward  manner  to  the  following: 
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pKyv|i(— H^)  = 


P«,v,i(-H')  = 


Pyyv|i(“H^)  “ 


*“  ■  «  Av,— n  ^  A— V,— 

“i^Bv.ji  +  B-v.n  •“  “■  B-v,-|i^ 

—i^Dv^  —  D-v,n  +  Dv,-|i  — 

-(Ev.,  -  E-v.,  -  Ev.-,  +  E-v.-,)]e-‘'-«' 


“  ^Av.g  +  A-v,j»  +  Av,-|i  +  A-v.-h^ 

-i(Dy.,  -  D.V.M  +  Dv.-,  -  D-v.-^) 

"“^Ev.|l  “  E— v.H  ”  Ev,— |l  E— V,— |l^jc 


+  V  kv(i  ^Av.fi  +  A-v,|4  +  Av,-|i  +  A-v.-fi^ 

+  ^  B— v,|l,  “  V,— 

+[Dv.^  -  D-».„  +  Dv,-g  -  D.v..,,) 

-i(Ev.,  -  E-V.K  -  Ev.-,  + 


“  '^'£2^  |l^^^Av,|i  ^  A-v,|i  +  Av»-^  +  A— V,— |i^ 

— i^Bv,|i  +  B— v,|i  —  Bv,-ii  *  B-.v,..|i^ 
”i^Dv,|i  **  D— v,n  +■  Dv.—ii  “  ^—v,— 

-[Ev.^  -  E-v.,,  -  Ev.-m  +  E-v.-M)]e-''-‘«' 


(40a) 


(40b) 


(40c) 


(40d) 
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+  +  A-v.n  +  Av.-i»  +  A-v.-n^ 

+  ^Bv,|i  ^  B— v.n  Bv,— |i  “  B— V,— 

+(Dv,„  -  D-v.,,  +  Dv.-m  -  D-v.-m) 
-i^Ev.„  -  E_v.|i  -  Ev.-n  +  E-v.-iijje"^’*”^ 


P«v,.(-HO  = 


•~k5„|^^Av.ii  +  A_»,ji  +  Av.-|»  +  A-v.-nj 
*i^Bv,|i  B— v,ii  ^  Bv.—ii  “  B— V,— 

-i^Dv.h  -  D-v.i»  +  Dv.-,i  -  D-v.-m  j 
-^Ev.(i  -  E_v.|i  -  Ev.-|»  +  E-v.-d^je"’'”'**^ 


(40c) 


Thus,  given  that  the  rectangular-harmonic  coefficients  have  previously  been  determined  as 
explained  in  section  3.2,  the  Fourier  coefficients,  which  are  associated  with  the  local 
rectangular-harmonic  potential,  are  ea^y  conq>uted  through  eqs.  (40a  through  40f)  and  that 
insoted  into  eqs.  (38a  through  38f),  which  may  be  used  in  conjunction  with  an  FFT  algorithm  to 
genoate  all  of  the  local  magnetic-gradient  conq>onents  of  the  surveyed  area.  Figures  4a  through 
4f  illustrate  the  result  of  this  con:q)utational  process  for  6  local  magnetic-gradient  componoits  of 
die  Northern  Juan  de  Fuca/Explora^  plate  region.  Figures  4a  through  4f  corre^nd  to  the 
magnetic  gradients  ,  BB^  ,  3B„  ,  BB^^ ,  BB^^ ,  and  BB^  ,  respectively.  Figure  4f  in  particular 

clearly  di^lays  die  nordi-south  oriented,  magnetic-reversal  lineation  patterns  of  the  Juan  de 
Fuca  ridge  axis  and  die  magnetically  quiet  subduction  zone  adjacent  to  the  axis. 
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de  Fuca  XX-Componmt  Magnetic-Field  Gradient;  Units;  nT/km 


Compon^mt  Magnetic-Field  Gradient;  Units:  nT/km 


Figure  4c.  Juan  de  Fuca  XZ-Componcnt  Magnetic-Field  Gradient;  Units;  nT/km 


4.  INVERSE  MODELING  (MAGNETIC  DEPTH-TO-SOURCE) 

The  previously  derived  rectangular-harmonic  model  of  the  northern  Juan  de  Fuca/Explorer 
plale  r^tm  is  of  phenomenological  origin  since  it  is  based  on  "observed"  survey  data. 
Consequently,  the  uniformly  ^aced  magnetic  component  grids  of  die  mi^netic-field  vector 
B^(x,y,z),  d^icted  in  Figures  3a  through  3c,  and  the  uniformly  spaced  magnetic  con^noit 
grids  of  the  magnetic-gradient  tensor  A^(x,y,z),  depicted  in  Figures  4a  through  4f,  all  of  which 

were  generated  from  the  model,  are  also  of  phenomenological  origin  and  are  dioefore  annotated 
with  the  subscript  "obs".  A  single  grid  location,  then,  has  attached  to  it  three  vector-magnetic 
conqKMient  observations  and  nine  (five  of  vdiich  are  indepoident)  gradient-tensor  component 
observations.  We  now  a^:  What  is  the  location,  distribution,  and  equivalent  crustal 
magnetization  M(r)  that  generates  the  observed  magnetic-vector  and  magnetic  gradient-tensor 
fidds? 

The  answo*  to  this  question  draws  upon  Maxwell's  equations,  which,  for  the  magnetostatic 
conditions  involved  in  this  problem,  reduce  to  the  following  Poisson  equation,  written  in 
Gaussian  units,  for  the  scalar  magnetic  potential  (^r); 

V»*(r)=4«V*M  (41) 

which  has  the  wdl-known  solution; 

V' 


udiere  V(r')  is  the  volume  occupied  by  the  magnetized  crustal  materials  and  whou  r  and  r' are 
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tile  poaticMi  vectors  from  the  origin  of  the  rectangular  coordinate  system  to  the  point  of 
observation  r(x,y,z)  and  to  the  location  r’(x',y’^')  of  the  magnetized,  oceanic  crustal  material 
within  the  volume  V  tiiat  generates  the  observed  fields.  Throi^  the  application  of  the 
divergence  theorem,  as  well  as  through  the  use  of  various  transformations  and  identities,  all  of 
\riiich  are  discussed  in  detail  by  Quinn  and  Shiel  (1993),  eq.  (42)  can  be  put  into  the  following, 
more  useful,  form: 

*(r)-  -V.  (43) 

V! 

Two  assunqitimis  are  now  made.  First,  the  magnetic-field  vector  and  gradient-tensor 
conqionents,  measured  at  the  point  r(x,y,z),  which  corresponds  to  a  particular  grid  location,  are 
generated  by  magnetized  sources  in  that  part  of  the  oceanic  crust  which  is  centered  directly 
below  the  grid  location  under  consideration.  Second,  the  distribution  and  strengtii  of 
magnetizaticm  in  the  oceanic  crust  is  slowly  varying  in  the  lateral  directions  away  from  this  grid 
location,  so  that  within  the  confines  of  a  few  kilometers  (e.g.,  S  km)  of  this  point,  the 
magnetization  may  be  considered  uniform  witiiin  the  volume  V*  and  hence  independent  of  r'. 
The  justification  for  the  first  assumption  rests  on  the  knowledge  that  the  vector  magnetic-field  in 
the  far-field  approximation  varies  as  R while  the  magnetic  gradient-tensor  field  varies  as  R~*, 
where: 

R  ■  I  r  -  r'  I  =  y(x  -  x')*  +  (y  -  +  (z  -  (44) 

is  tiie  distance  between  a  point  within  the  magnetized  source  volume  V  and  the  point  of 
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observation.  The  second  assunqition  rests  on  the  knowledge  that  the  peak>to-pcak  sepaiatiKMi  of 
magnetic  reversals,  at  least  in  the  Juan  de  Fuca  region,  is  on  the  order  of  60  km.  Consequently, 
eq.  (43)  may  be  written  as: 


4Kr)  =  -M.vJ  j^dV  (45) 

V' 


Taking  the  negative  gradient  of  this  potential  yidds  the  magnetic-field  vector 


B(r).M.vJ  (4«) 

yt 


The  volume  V  may  be  characterized  as  a  single  rectangular  prism  of  dimoisions  x  Xy  x  , 
v^ch  were  a  priori  taken  to  be  S-km  x  S-km  x  l.S-km  on  the  basis  of  drill  hole  data  and  our 
knowledge  of  die  far-field  bduvior  of  the  magnetic  field  and  its  gradients.  The  prism’s  coiter  is 
located  directly  below  the  known  observation  point  (Xq  ,  y^),  at  an  unknown  source  depth  Zg , 
vtdiich  must  be  determined  from  the  observed  data  along  with  the  magnetization  vector  M,  which 
has  constant  conqxments  M,  =  M, ,  Mj  =  M, ,  and  M,  =  M, . 

Using  the  rectangular-prism  geometry,  eq.  (46)  takes  the  following  form  in  tensor  notation: 

B,i(x,y,2)  *  AjuM*-  (47) 

whoe  Einstein  summation  over  the  repeated  index  X  (one  raised  and  one  lowered)  is  assumed 
rad  y/hat  the  elements  of  the  A  matrix  are  given  by  the  relation: 
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A»;i(x.y,z) 


L-^ 

J 

Jy.-- 

UJ' 

(48) 


These  mstrix  dements  can  be  evaluated  analytically  as  discussed  by  Quinn  and  Shid  (1993). 
The  results  of  diis  evduadon  are  given  in  Table  5,  whoe  the  notation  |.< ,  y> ,  ^  indicates  that  the 
parameters  x' ,  y* ,  and  z*  are  to  be  evaluated  at  dieir  req)ective  iq>per  and  lower  limits.  That  is: 


I,'. 


y'a' 


iy«* 


.it 


(49) 


Consequently,  each  matrix  elemoit  in  Table  S  has  dght  terms.  Furdiermore,  it  is  clear  from 
eq.  (48)  that  die  A  matrix  is  symmetric.  It  can  also  be  diown  that  this  matrix  has  zero  trace  if 
die  point  of  observadmi  is  outside  the  volume  V  and  a  trace  equal  to  4x  if  the  point  of 
observadon  is  inside  that  volume.  This  matrix  is  also  a  tensor. 

Taking  the  gradient  of  eq.  (47)  yields  the  magnetic-gradient  tensor  due  to  a  uniformly 
magnetized  prism: 


^)iv(x,y,z)  =  Agx/vM^  (50) 

\^iere  the  derivatives  of  the  A  matrix  elonents  can  also  be  evaluated  analytically  and  are  givoi 
in  Table  6. 

The  nu^etization  of  die  prism  can  be  formally  determined  by  inverting  eq.  (47),  vdiich 
yields  die  result 
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Tables.  Elonoitsof  die  A  Matrix 


Table  6.  Derivativcf  of  the  A  Matrix  Eloaents 

An/i  *  +  t[tS^  + 

An/i  =  " 

An/j  =  - 

An/I  »  ”  t[i7^o] 

An/2  *  -  r[iT^o]  ^ 

An/3  *  ~  “r  li'.y'.i' 

An/I  =  “ 

An/2  =  “  ■r  I  »'.y'.*' 

An/3  =  -  l’[ii+(y^] 

A21/1  =  -  7[i7^]  l»'.y'.*' 

Aii/2  *  -t[i7^] 

A21/3  *  “  X  ln'.y'.i' 
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Table  6.  Derivatives  of  the  A  Matrix  Elemeuts  (con.) 


I 

I 


Table  6.  Derivatives  of  the  A  Matrix  Elements  («>n.) 


wdiere  the  matrix  X  is  the  inverse  of  A .  Then,  inserting  this  result  into  eq.  (SO),  yields: 


This  is  the  main  result.  It  is  a  nonlinear  equation  for  the  magnetic-source-depth  Zq,  givoi  the 
"observed"  magnetic-field  vector  and  gradient-tensor  components  at  the  point  (x^  ,  y^.  The 
nu^etic-source-depth  can  be  detenmned  by  a  nonlinear  least-squares  mininuzation  of  die 
chi-square  )  function: 


X*-  —  obt)^ 

|ta|  vsl 


with  reflect  to  Zg  .  For  the  Juan  de  Fuca  region,  this  chi-square  fimction  was  evaluated  as  the 
value  of  was  varied  at  0. 1-km  intervals  fiom  the  known  bathymetric  surface  to  slightly  beyond 
the  estimated  Curie  depth,  which  in  turn  was  based  on  a  spherical-harmonic  model  of  heat 
flow  (Pollack  and  Quqmian  [1987]).  The  magnetic-source-depth  then  correqionds  to  the 
minimnm  value  of  .  Performing  this  procedure  at  each  grid  point  of  die  survey  area  yields  the 
magnedc-source-dqith  of  the  region  as  depicted  in  Figure  S.  Clearly  defined  in  this  figure,  to  a 
hi^  depee  of  resolution,  is  the  subduction  zone  just  west  of  Vancouver  Island  and  west  of  the 
state  of  Washington. 
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Having  evaluated  the  magnetic-source-depth  at  each  grid  point,  the  A  matrix  and  its  inverse 
are,  as  a  consequence,  also  conqiletely  determined  at  each  grid  point.  This  allows  the  prism 
magnetization  to  be  evaluated  at  each  ^d  point  through  the  linear  inversion  described  by 
eq.  (SI).  The  resulting  three  conqKments  of  magnetization,  which  in  the  Gaussian  units  used 
here  are  expressed  in  lumoTesias,  at  each  grid  point  of  the  Northern  Juan  de  Fuca/Explorer  plate 
survey  region  are  illustrated  in  Figures  6a,  6b,  and  6c.  Figures  7a  through  7e  are  selected 
magnetic-source-depth  profiles  of  the  Z-corxq>onent  magnetization.  Figures  7a  and  7b  are 
profiles  along  constant  longitude  bands,  while  Figures  7c  through  7e  are  profiles  along  constant 
latitude  bands.  The  depths  and  magnetizations  in  these  figures,  though  reasonable,  are  not 
considered  to  be  absolute,  since  the  a  priori  assuii^)tions  are  only  approximately  true.  Secondly, 
in  order  to  speed  up  the  inversion  process  an  extra  "nonessential"  constraint  was  applied  to  the 
total  magnetization,  which  limited  its  value  to  4  Amps/meter  in  SI  units.  The  total  magnetization 
conqiuted  during  the  inversion  always  reached  diis  limit,  indicating  that  our  choice  was  too  low. 
In  the  mid  to  late  1980's,  when  this  work  was  actually  performed,  this  limit  was  thought  to  be 
reasonable.  More  recent  dredge  sampling  indicates  that  the  magnetization  of  younger  basalt 
materials  near  ridge  crests  is  as  much  as  an  order  of  magnitude  greater.  The  vector  components 
were  allowed  to  vary  freely  within  the  limits  imposed  by  the  constraint  on  the  total 
magnetization.  The  relative  depths  and  magnetizations,  so  derived,  appear  at  least  qualitatively 
correct  but  may  require  appropriate  scaling  to  be  quantitatively  useful.  The  fluctuations  observed 
in  the  depth-to-source  magnetization  profiles  are  generally  consistent  with  serpentinization  faults 
goioally  found  at  ridge  crests  (Francis,  1981).  For  a  variety  of  reasons,  we  have  not  found  the 
resources  to  rerun  these  computations  without  this  extra,  nonessential  constraint  on  the  total 
magnetization  as  we  would  prefer  to  do. 
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Figure  6b.  Juan  de  Fuca  Y-Component  Magnetization;  Units:  nT 
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Figure  7a.  Juan  de  Fuca  Z-Component  Magnetization  Profile  (nT)  Along  Longitude  129.65W 


Figure  7b.  Juan  de  Fuca  Z-Component  Magnetization  Profile  (nT)  Along  Longitude  129. SOW 
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Figure  7c.  Juan  de  Fuca  Z-Component  Magnetization  Profile  (nT)  Along  Latitude  47.30N 


Figure  7d.  Juan  de  Fuca  Z-Component  Magnetization  Profile  (nT)  Along  Latitude  48.23N 


Figure  7e.  Juan  de  Fuca  Z-Component  Magnetization  Profile  (nT)  Along  Latitude  48.40] 


5.  COMMENTS  ON  NONUNIQUENESS 


The  invetse  gec^otential  modeling  problem  is  notoriously  ntmunique.  That  is,  for  a  given  set 
of  nu^snetic-fidd  survey  measurenaents,  many  disnibutions  of  magnetization  may  be  found  diat 
can  rqproduce  die  observed  magnetic  field.  This  is  true  regardless  of  the  accuracy  of  the 
mi^netic-field  measurements.  The  way  around  this  difficulty  is  to  make  reasonable  assunqitions 
based  tm  a  priori  infbrmatitm  that  may  restrict  die  possible  source  geometry  and  reduce  die 
mmiba*  of  possible  udmown  parameters.  In  die  present  case,  we  assumed  a  rectangular-imsm 
geometry  based  on  the  idea  dut  die  lateral  magnetic  character  of  a  large  oceanic  region  does  not, 
as  a  rule,  diange  drastically  over  die  course  of  a  few  kilometers.  Furdtermore,  this  choice  of 
geometry  was  based  on  drill  hole  and  other  data  which  indicated  that  the  largest  source  of 
magnetization  in  the  ocean  crust  is  confined  to  the  basalt  laya.  The  choice  of  the  rectangular- 
prism  geometry  still  left  nine  parameters  to  determine,  three  conqionents  of  magnetization,  three 
coordinates  qiecifying  the  center  location  of  the  prism,  and  diree  ^ism  dimension  parameters. 
Then,  given  tmly  the  vector  magnetic-field  measuronents  (i.e.,  no  gradient  data)  as  is  frequoidy 
the  case,  defining  the  elements  Cj  ,  j  =  1,2,...6  to  corre^nd  to  the  prism  parameters 

and  \ )  and  neglecting  magnetic-field  measurement  errors,  we  find  that  a  variation  of  the 
parameters  in  eq.  (47)  yields: 


5B^  *  '  M*'6cj  +  A|ix5M^  =  0 


(54) 


which  inqilies: 

•Cj  ^ 


(55) 
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Eqiiations  (54)  and  (SS)  state  diat  without  any  further  constraints  on  the  geometric  parameters 
Cj ,  any  variation  in  one  or  more  of  these  parameters  can  be  offset  by  a  compensating  variation  in 
the  magnetization,  ^^ch  will  leave  the  magnetic  field  unchanged.  So,  without  fiirtha 
constraints,  the  solution  is  ambiguous.  This  ambiguity  is  eliminated  by  first  demanding  that  the 
prism  be  centered  directly  below  the  particular  grid  point  under  consideration.  This  means  that 
5c,  =  6c2  =  0  .  Furthor  noting  that  the  magnetic  field  falls  off  sharply  with  distance  from  the 
source,  a  cutoff  length  to  the  lateral  dimensions  of  the  prism  can  be  imposed  (we  used  5  km),  so 
that  Sc^  -  Scj  =  0.  Drill  hole,  seismic,  and  other  data  indicate  that  in  the  ocean,  the  primary 
source  of  the  crustal-generated  magnetic  field  is  located  in  the  basalt  layer  which  is 
q>proximately  l.S-km  thick.  Using  this  value  for  c^  sets  Sc^  -  0.  This  leaves  the  source  depth 
c,  -  2^ ,  as  the  remaining  undetermined  geometric  parameter,  which,  without  further  constraining 
information,  will  still  lead  to  a  nonunique  solution.  The  required  new  information  comes  from 

our  knowledge  of  the  magnetic  gradient.  A  variation  of  the  geometric  parameters  in  eq.  (52), 

again  ignoring  measurement  errors,  so  that  5^^^  =  0,  yields: 

+  AHX/v2g^B,}5cj  =  0  (56) 

Since  5<^  =  0  for  j  ^  3  ,  eq.  (56)  reduces  to: 

=  0  (57) 


This  is  actually  a  set  of  9  equations,  5  of  which  are  indepoident.  Each  of  these  equations 


implifta  that  the  source  depth  uncertainty  SC3  =  Sz^  =  0,  since  the  term  within  the  brackets  is  not 
zero  in  general. 

Thus,  the  mediod  used  to  determine  the  magnetizations  and  source  depths  yields  a  unique 
solution  at  least  to  the  extent  that  the  arbitrarily  inqMsed  constraints  on  the  lateral  prism 
dimensitms  are  satisfied.  The  largest  source  of  error  in  this  procedure  results  from  the 
assunqrtion  that  just  one  magnetic  source  exists  beneath  a  q)ecified  grid  location.  Actually, 
weaker  sources  exist  at  shallowo*  depths.  The  characteristics  of  these  weaker  sources  can  be 
determined  through  an  iteration  process  which  removes  the  magnetic  fields  generated  by 
ma^etizations  obtained  through  the  inversion  procedure  from  the  original  survey-generated 
rectangular  grids  of  crustal  magnetic-field  and  gradient-toisor  residuals.  This  process  may  be 
repeated  until  the  residual  magnetic-field  vector  is  of  the  order  of  the  measurement  errors 
associated  with  the  survey.  For  each  iteration,  the  same  uniqueness  arguments  apply,  as  do  the 
provisos. 

As  a  final  note,  it  should  be  kept  in  mind  diat  we  could  have  minimized  eq.  (53)  with  respect 
to  die  latoal  dimensions  of  the  prism  (X,  and  X^)  as  well  as  the  prism's  source  depth  Zg.  This 
removes  the  arbitrary  nature  of  the  lateral  dimensions  at  the  expense  of  considerable  amounts  of 
conqiuter  time.  However,  this  would  not  resolve  the  ambiguity  problem  since,  as  the  above 
exercise  indicates,  for  any  particular  source  depth,  an  infinite  set  of  lateral  dimension 
combinations  will  satisfy  the  minimization  problem.  The  addition  of  one  more  physically 
reasonable  constraint  would  resolve  the  ambiguity.  For  exanqile,  one  could  also  simultaneously 
damand  that  the  prism  volume  be  a  minimum.  The  physical  basis  for  this  particular  constraint  is 
radier  dubious  and  no  other  physically  reasonable  constraints  were  found.  So,  this  approach  was 
not  taken. 
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